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Abstract 

The presence of charge promotes the propensity of a rotating black hole to support the prop- 
agation of electromagnetic plane waves with negative phase velocity (NPV) in its ergosphere, 
whether the Kerr-Newman or the Kerr-Sen metric descriptions of spacetime are considered. 
Striking differences in the NPV characteristics for the Kerr-Newman and Kerr-Sen metrics 
emerge from numerical studies, particularly close to the outer event horizon when the magni- 
tude of the charge is large. 
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1 Introduction 

Electromagnetic plane waves with negative phase velocity (NPV) have been the centre of much 
attention lately [1]. They give rise to many exotic phenomenons — the most prominent being 
negative refraction [2] — which are not associated with their positive phase velocity (PPV) coun- 
terparts. Much of the interest in NPV propagation has stemmed from the realization of negatively 
refracting metamaterials [3] . In addition, the prospects of NPV propagation in special [4] and gen- 
eral [5] relativistic scenarios, and the concomitant implications for observational astronomy, have 
prompted study and speculation within the past few years [6] . 

It is now well established theoretically that, for certain spacetime metrics, NPV propagation 
can be supported in free space [7, 8]. In particular, NPV can arise in the ergosphere of a rotating 
black hole described by the Kerr metric [9] . The NPV-supporting regions are concentrated towards 
the equator of the ergosphere, and this concentration becomes exaggerated as the angular velocity 
of the black hole is increased. Here we address the question: how does the presence of charge effect 
the propensity of a rotating black hole to support NPV propagation? 

The most commonly encountered spacetime description of a charged rotating black hole is 
provided by the Kerr-Newman metric [10]. This solution to the Einstein field equations was derived 
in the 1960s [11, 12]. In 1992, Sen developed an alternative metric solution describing a charged 
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rotating black hole, based on the action of a heterotic string [13]. Since then, various aspects of the 
structure and properties of this Kerr-Sen metric have been established [14, 15, 16]. In the following 
sections, we consider both the Kerr-Newman and the Kerr-Sen metrics in our investigation of the 
NPV characteristics of a charged rotating black hole. 



2 Theoretical framework 
2.1 Preliminaries 

Electromagnetic propagation is considered in free space. The spacetime curvature is described by 
the metric g a ^ with signature (+,—,—,—). Here we outline the theory of gravitationally assisted 
NPV propagation, comprehensive details being available elsewhere [5, 9]. 

Following the approach originally put forth by Tamm [17, 18, 19, 20], we exploit the fact that 
electromagnetic propagation in curved spacetime in free space is equivalent to propagation in a 
fictitious, instantaneously responding, bianisotropic medium in flat spacetime, characterized by the 
constitutive relations 

Di = je m E m + El mn V m H n I 

r • (1) 

B(, — JimHm — ££ mn T m E n J 

Herein, Eg, B^, Dp, and are the components of the conventional electromagnetic field vectors; 
£imn is the three-dimensional Levi-Civita symbol; and the components of the tensor j£ m and the 
vector T m are defined as 



I \l/2 9 

r _ 90m 

m — 

5oo 

In order to consider planewave propagation, we rewrite (1) in 3x3 dyadic/3 vector form as 

D(ct,r) = e 7(ct,r) • E(ct,r) - -T{ct,r) x H(ct,r) 

B(ct,r) = Hol(ct,r) ' H(ct,r) + -Tfct.r) x E(ct,r) 
= c 



(2) 



(3) 



wherein j(ct,r) is the dyadic-equivalent of 7^ m , T_(ct,r) is the vector-equivalent of T m , the scalar 
constants e and fj, denote the permittivity and permeability of vacuum in the absence of a gravi- 
tational field, and SI units are adopted. In source-free regions, E_, B_, D, and H_ are related by the 
Maxwell curl postulates as 

d 

V x E(ct,r) + ^B(ct,r) = 



d 

V x H(ct,r) - —D(ct,r) = 



(4) 



Let us emphasize that, while the spatial (r) and temporal (t) coordinates have been separated in 
(3) and (4), the underlying spacetime remains curved. 

Our aim is to explore electromagnetic planewave propagation in the spacetime of a black hole 
of geometric mass rribh and charge qbh, rotating with angular momentum a^h about the Cartesian z 

§ Greek indexes take the values 0, 1, 2 and 3, Roman indexes take the values 1, 2 and 3, x° = ct, where c is the 
speed of light in vacuum in the absence of all gravitational fields, and a: 1 ' 2 ' 3 represent the three spatial coordinates. 
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axis, in conventional units. Two different descriptions of this spacetime are considered, as provided 
by the Kerr-Newman metric and the Kerr-Sen metric. The line elements for these are convention- 
ally expressed in terms of Boyer-Lindquist coordinates. In order to establish the corresponding 
forms of the dyadic 7 and vector T, the Kerr-Schild Cartesian coordinate representations for these 
metrics are presented in the following two subsections. 



2.2 Kerr— Newman spacetime 

The Kerr-Newman line element is given in Boyer-Lindquist coordinates as [21, 22] 



ds 2 = 



A 



KN 



R 2 + a bb cos 2 
R 2 + a 2 bb cos 2 9 

where R is given implicitly via 



[dt — a bb sin 2 Odcf)' 



sin 



R 2 + al h cos 2 6 



[(R 2 + a 2 bh ) d(j)-a hb dtf 



dR 2 - (R 2 + a 2 bb cos 2 9) d9 2 



R 2 



2 1 2 , 2 
x +y + z 



-(5) 



and 



A 



KN 



R 2 - 2m bb R + a 2 h + q 2 b 



"bh 1 Hbti ■ 

By applying the coordinate transformation (t, (f>, R, 9) — > (t, x, y, z) specified by 

(i? cos 4> + a bb sin cf) sin 9 , 
(i? sin cj) — a bb cos (f) sin 9 , 
R cos 9 , 

dt = dt + dR- 

^bh 



X 

V 

z 



R2 + < <t R, 



d(j) = d(f> 



A 



-dR, 



KN 



(5) 

(6) 

(7) 

(8) 
(9) 
(10) 

(11) 
(12) 



the line element (5) is expressed in Kerr-Schild Cartesian form as 
df - dx 2 - dy 2 - dz 2 - {2mhbR " ql) ^ * 
R 



ds 2 



R 4 + 



fluu Z 



dt — —dz 

R R 2 + a 2 h 



(xdx + ydy) 



(xdy — ydx) 



(13) 



R 2 + a 2 bb 

The corresponding metric components g a p — from which the associated dyadic 7 and vector T are 
immediately delivered by the definitions (2) — are provided explicitly in Appendix 1. 

We focus on the regime m bh > a bb . Furthermore, our interest lies in the ergosphere demarcated 
as R + < R < Rs+! with R = R + representing the outer event horizon and R = Rs+ representing 
the stationary limit surface [23]. The outer event horizon occurs where Akn = 0. Thus, 



R + =m bb+ ^ 



<?bh- 



The stationary limit surface R = Rs + occurs where goo = 0. Thus, 

Rs+ = m bh + Mm 



bh 



9 ( a bh Z \ 



(14) 



(15) 
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2.3 Kerr-Sen spacetime 

We now turn to the Kerr-Sen line element which is given in Boyer-Lindquist coordinates as [14, 16] 



ds 2 



S 2 + a 2 cos 2 6 



(dt — a bh sin 2 Odcf) ' 



sin 



S 2 + a 2 cos 2 9 



(s 2 + a 2 h ) d(f> - a bh dt 



S 2 + a 2 h cos 2 9 



-dS- 



-( 



where 



S 2 + a£ h cos 2 )de 2 , 



P 



A KS = S'-2 (^m bh - - ) S + a z bh , 

S 2 = s 2 + ps, 



= 



and S satisfies 
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52 2,2,2 2 

<S = x + y z + z z - a b 
Under the coordinate transformation (t, 4>, S, 6) — ► (t, x, y, z) specified by 

S cos (f> + a bh sin cf) ) sin , 



z 
dt 
d4> 



[s sin (j) — a bh cos cj>J sin ( 
S cos 9 , 



= dt + dS- 



S 2 + ^ 



^bh 



the Kerr-Schild Cartesian form of the line element (16) emerges as 



ds' 



dt - dx 2 - dy 2 - dz 2 + 



(3 2 (s 2 + a 2 h cos 2 6) ^ _ 2mhJ 2 S 
(jP + 4S 2 ) A KS S 4 + < z2 



dt — -dz — . 

S S 2 + a 2 bh 



(xdx + ydy) - ~ ° bh (xa?y - ydx) 



S 2 + a 2 



bh 



(16) 

(17) 
(18) 
(19) 

(20) 

(21) 
(22) 
(23) 

(24) 
(25) 



(26) 



In Appendix 2, explicit expressions for the corresponding metric components g a/ 3 are presented. 
From these, the associated dyadic 7 and vector T immediately delivered by the definitions (2). 

In analogous manner to the case for the Kerr-Newman metric, the outer event horizon arises 
at S = S+ with 



6+ = m bh - - + 



m bh 



(27) 
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whereas the stationary limit surface arises at S = Ss + with 



i>s+ = m bh - - + 



ft 



Ss 4 



(28) 



As in §2.2, the ergosphere is the region bounded by the outer event horizon and the stationary limit 
surface; i.e., S+ < S < Ss + - 

2.4 Planewave propagation 

Now we turn to the propagation of plane waves in the medium characterized by (3). By choosing a 
sufficiently small neighbourhood 1Z specified by the Cartesian coordinates (x, y, z), we approximate 
the non-uniform metric g a p with the uniform metric g a p throughout 1Z [7]. Within the neighbour- 
hood 1Z, we utilize the uniform 3x3 dyadic 7=7 and the uniform 3 vector T = Tj^, where 

7 and T are calculated from (2). Due to the uniform approximation, attention is restricted to 
electromagnetic wavelengths that are small compared to the linear dimensions of 1Z. 
The three-dimensional Fourier transforms 



/•oo /»oo roo 

E(ct,r) = - / / / E(cj/c, k) exp [i(k • r — tut)] dw dk\ dk,2 
c J— 00 J— 00 J— 00 

Y poo rco roo 

H(ct,r) = - / / / / c, k) exp [i(k • r — ujt)] dco dki dk2 

c J —00 J —00 J —00 



(29) 



decompose the electromagnetic fields within 1Z, Here, the wavevector k is the Fourier variable 
corresponding to r, i = \f—l and uo is the usual temporal frequency. By combining (4) and (29), 
and then solving the resulting 4x4 eigenvalue problem, the wavevector component k^ is determined 
[24]. Only propagating (non-evanescent) plane waves are considered; i.e., k% 6 M. 

NPV propagation occurs when the phase velocity vector casts a negative projection on to the 
time-averaged Poynting vector — i.e., 



k ' (P(uj/c,k)) t <0. 



(30) 



where the time-averaged Poynting vector (P_(uj/c, k) )t is calculated from the complex-valued pha- 
sors of the electric and magnetic fields, E(uo/c, k) and ^i.(i^/c,k). The left side of (30) is given by 
[9] " " 



k'{P(u/c,k)) t = 2^J^j (l^ Q ! 2 g« ' 1 ' ga + I A| 2 gb 'l'2b) k'g'p. 



(31) 



The complex-valued amplitudes A a ^(co/c, k) are fixed by initial and boundary conditions, whereas 
the unit vectors 



7 1 • w 

|7 _1 • w\ 

-1 . 



> . 



7 



P X §aj 



It 1 * (p x § a ) I j 



(32) 
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For the numerical results presented in the following section, we have chosen 



w = < 



2^5 + (Px+PyY 



(Pz,Pz, ~Px ~ Py) for p z / 



(33) 



+ p2 



(Py,~Px,0) 



for p z = 



with p = (p x ,Py,Pz) m the Cartesian coordinate system. 

For both the Kerr-Newman and Kerr-Sen metrics | 7 | > 0. Hence, the sufficient condition for 
NPV propagation 

P a <0 

n<o 



(34) 



arises, wherein 



Pa = (§ a «2 [k'X'P 

Pb = (§ 6 'g *e 6 , ) (fc *2 *e) 

The two wavenumbers k = k^ given by 



(35) 



± 



a; 



A; 



■ 7 • f ± ^(i'7'f) 2 -fc'7'fe^f'7'f-|7|^ 



V 



• 7 • 

develop for the arbitrarily oriented wavevector k = kk with k = (sin 9 cos 0, sin 6* sin 4>, cos 0). 



(36) 



3 Numerical results 

The NPV characteristics of an uncharged rotating black hole, as described by the Kerr metric, 
have been established previously [9]. In this section, we report on the effect that charge has on the 
propensity for a rotating black hole to support NPV propagation in its ergosphere. This study was 
carried out by numerically evaluating (31), for both the Kerr-Newman and Kerr-Sen descriptions 
of the blackhole spacetime. 

The following strategy was implemented: for fixed values of m bh , a hh and q bh , we considered 
surfaces of constant R for the Kerr-Newman metric, and surfaces of constant S for the Kerr-Sen 
metric, with all surfaces being contained within the respective blackhole ergospheres. The prospects 
for NPV propagation for all propagation directions k = (sin cos </>, sin 6 sin 0, cos 9) were calculated 
over a grid of points on surfaces of constant R and S. At each coordinate point (x, y, z) on the 
grid, the quantities P a and Pb were computed at 1° increments for 9 G [0°, 180°) and (f> G [0°,360°), 
for both k = k + and k = k~ . For each value of (9,(p) at which NPV was found to be supported 
(i.e., P a < and P5 < 0), a score of one point was assigned to the (x,y,z) grid point. The total 
score at each grid point was translated to a colour for the convenient presentation of results. 

Figure 1 shows the NPV propagation maps for Kerr-Newman spacetime for a charged rotating 
black hole with a bh = m bt ^3/4. Three constant R surfaces within the ergosphere are considered, 
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namely R = 0.99i?+ + OMR s+ (in Figures 1(a) and (d)); R = 0.5R+ + 0.5R S+ (in Figures 1(b) 
and (e)); and R = 0.01-R+ + 0.99i?s + (in Figures 1(c) and (f)). For each surface, a small value 
of q bb (Qbh = 0.1 in Figure 1(a), (b), and (c)) and a large value of q bh (q bh = 0.45 in Figure 1(d), 
(e), and (f)) are chosen. We note that q bh < 0.5 in order for the ergosphere boundaries to remain 
real- valued. 

In Figure 1(a)— (c) that, for q hh = 0.1, the regions of NPV propagation are most heavily con- 
centrated near the equator of the ergosphere and NPV becomes less likely towards the poles. Also, 
NPV propagation is more likely near the outer event horizon and less likely near the stationary 
limit surface. Increasing the value of q hh does not markedly change the qualitative nature of the 
NPV characteristics, as illustrated by Figure l(d)-(f). However, the incidence of NPV propagation 
is considerably greater at q bh = 0.45 as compared with q bh = 0.1, at all three constant R surfaces. 

The corresponding NPV propagation maps calculated using the Kerr-Sen metric are presented 
in Figure 2. A comparison of Figure 2(a)-(c) with Figure l(a)-(c) reveals that the NPV charac- 
teristics of the Kerr-Sen and Kerr-Newman metrics are quite similar at q bb = 0.1. However, stark 
differences between the Kerr-Sen and Kerr-Newman NPV characteristics emerge at q bh = 0.45, 
particularly near to the outer event horizon. For the Kerr-Sen metric, we see in Figure 2(d) that 
the distribution of NPV propagation is concentrated in a region between the ergosphere pole and 
equator, while a relatively low incidence occurs at the equator. This is quite different from the 
situation for the Kerr-Newman metric that is displayed in Figure 1(d). 

We note that the NPV propagation maps at q bh = 0.1, for both the Kerr-Newman metric and 
Kerr-Sen metric, are quite similar to those reported for the uncharged rotating black hole [9]. This 
reflects the fact that both the Kerr-Newman metric and the Kerr-Sen metric become identical to 
the Kerr metric in the limit q bh — ► 0. Furthermore, regardless of charge, NPV is not supported at 
the ergosphere pole. 

4 Concluding remarks 

Our investigation has revealed that the propensity for a rotating black hole to support NPV prop- 
agation is highly sensitive to the presence of charge in the black hole. Specifically, 

• the overall incidence of NPV increases as q bh increases, for both the Kerr-Newman and Kerr- 
Sen metrics; 

• as q bh increases, differences in the NPV characteristics associated with the Kerr-Newman and 
Kerr-Sen metrics emerge; and 

• the greatest differences between the NPV characteristics of the Kerr-Newman and the Kerr- 
Sen metrics arise close to the outer event horizons of their respective ergospheres. 

These findings suggest the possibility of a NPV-based experiment to distinguish between the Kerr- 
Newman and Kerr-Sen metric descriptions of a charged rotating black hole. 
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Appendix 1 

The components of the Kerr-Newman metric obtained from the line element (13) expressed in 
Kerr-Schild Cartesian form are 

5oo = 1 - $KNi (37) 
5oi = + 5 K n)R 3 xt 2 + A KN {-a hh 5 K N[R 3 (a hh x + Ry) + a hh (-Rx + a hh y)z 2 } + R 3 xt}), (38) 

502 = cr((-l + 5 K N)R 3 yT 2 + A KN {a hh 5 K N[R 3 (Rx - a bh y) + a bh z 2 (a bh x + Ry)] + R 3 yr}), (39) 
503 = ^-{A KN R 4 + a 2 bh A KN [-(-l + 5 KN )R 2 + 8 KN z 2 ] + (-1 + 5 kn )R 2 t 2 }, (40) 

5n = -<?- 2 ((-l + 5 K n)R 6 x 2 t 4 - 2A KN R 3 xT 2 {a bh 5 K N[R 3 (a bh x + Ry) + a bh z 2 (-Rx + a hh y)\ 
-R 3 xt} + A 2 KN {al h 5 K N[R 3 (a bh x + Ry) + a bh z 2 (-Rx + a bh y)] 2 

+ (R 8 - R 6 x 2 + 2a 2 bh R 4 z 2 + a 4 bh z 4 )r 2 }), (41) 

2 ,R % xy{AKN — t) 2 rri n . . 9i 

512 = o- ( ^ \-6 K N{a hh A KN [R (a bh x + Ry) + a hh (-Rx + a hh y)z ] 

-R 3 xT 2 }{a bh A KN [R 3 (Rx - a bh y) + a bh (a bh x + Ry)z 2 ] + R 3 yr 2 }), (42) 
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513 = ^^(R 5 x(A KN - t) 2 t 2 - 5 KN [a 2 bh A KN (-R 2 + z 2 ) + R 2 T 2 ]{-a bh A KN [R 3 (a bh x + Ry) 
+a hh {-Rx + a hh y)z 2 ] + R 3 xt 2 }), (43) 

522 = -a 2 (a 2 bb S KN A 2 KN [R 3 {Rx - a bh y) + a bh (a hh x + Ry)z 2 } 2 + A KN {R 6 [A KN (R - y)(R + y) 
+2a hh 5 K Ny(Rx - a hh y)] + 2al b R 3 [A KN R + 5 K Ny(a bh x + i?y)]z 2 }a 4 h A,ft-Ar2: 4 )T 2 
+2A KN R^y 2 T 3 + (-1 + 5 KN )R%\\ (44) 



523 = —(R b y{A KN -TYT z -5 KN [al h A KN {-R z + z z ) + R z T z } 

x{a bb A KN [R 3 (Rx - a hh y) + a bh (a bh x + i?y)z 2 ] + R 3 yr 2 }), (45) 



-1 

+2A^#V[# 4 + a 2 h (i? 2 - ^ii 2 + S kn z 2 )]t 2 + (-1 + fev)#Vr 4 }), (46) 



533 = -^(a 2 T 2 {A 2 KN (R-z)(R + z)[R 8 + ai h (-l + 8 KN )R 2 z 2 -ai h S KN z 4 



and 5 a/3 = gp a . Also, 



_ fi 2 (2m bhJ R - q 2 ) 
SKN ~ R* + a 2 bb z* ' (47) 



r = a 2 h + J R 2 , (48) 
°" = A^r^ + a 2 ^)- (49) 

Appendix 2 

The components of the Kerr-Sen metric obtained from the line element (26) expressed in Kerr- 
S child Cartesian form are 

500 = 1 - Sks, (50) 

5oi = ${(-1 + S KS )S 3 x[-2a 2 bh A KS S + v(-A KS + pi/)] + 2A^[-a bh fey(S 4 + a 2 bh z 2 ) 

+x(S 4 + a 2 hh 5 KS z 2 )S}}, (51) 

502 = ${(-1 + fe)S 3 y[-2a 2 h A^ + K-A^s/3 + H] + 2A^ s< S[a bh fez( < S 4 + a 2 bh z 2 ) 

+ y(S 4 + a 2 hh 8 KS z 2 )S}}, (52) 

503 = ^[A KS z{2{S A + a 2 h [-(-l + 5 KS )S 2 + 6 KS z 2 ]} - (-1 + 5 KS )S(a 2 bh + 5 2 )/3) 

+ (-l + fes)S(a 2 h + S 2 )zH, (53) 
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&{S*a?[-(-l + 5 KS ){2a 2 bb A KS S + v[A KS (3 - pv}} 2 + A 2 KS (2a 2 bh S + (3v) 2 Q 
+4A KS S[a bb 5 KS S 3 xy{S 4 + a 2 bb z 2 )[-2a 2 bh A KS S + v{-A KS f3 + 2Sv + (3u)] 
-[a 8 bh A KS z 4 + a 6 hh A KS z 2 (2S 4 + 5 KS y 2 z 2 ) + a 4 bh A KS S 4 [S 4 + 28 KS {-x 2 + y 2 )z 2 } 
+S 7 x 2 u[pu - A KS p(l + Q] + a 2 bh S 3 (5 KS x 2 z 2 pu 2 - A KS {5 KS x 2 z 2 pv 
+S 5 [S KS y 2 + 2x 2 (l + ()}})]§ + 2al h 5 KS A KS xyz 2 (S 4 + a 2 hh z 2 )S 2 - A KS [2a 2 hh S 8 
+2atz 4 + a 4 bh (AS 4 z 2 + 5 KS x 2 z 4 ) - S 8 x 2 (l + C)]S 3 - A KS (S 4 + a 2 bh z 2 ) 2 S 5 }} , (54) 

= A 2 KS ^ 2 f Xyp2 f 2 KS - v)2 + S«xyp 2 C + -±- 2 (6 KS (S 3 x{2a 2 hb A KS S 

^KS ^KS^ 

+v[A KS P - HI + 2a bh A KS S(y(S 4 + a 2 bb z 2 ) - a bb xz 2 S)){S 3 y{-2a 2 bb A KS S 
+u[-A KS (3 + pu}} + 2a bb A KS S[x(S 4 + a 2 bb z 2 ) + a bb yz 2 S}})}, (55) 

--^—[5ks{A K sS 3 (3 + a 2 h A^[-2z 2 + Sp] - S(a 2 hh + 5 2 )pz,}(5 3 x{2a b 2 1 A^S 
+v[A KS f3 - pu\} + 2a bh A KS S[y(S 4 + a 2 h z 2 ) - a bh xz 2 £])]}, (56) 



-$ 2 {S 6 y 2 {-A 2 KS (2a 2 bb S + (ivfC + ("1 + S KS )[2a 2 bh A KS S + v(A KS (5 - up)] 2 } 
+AA KS S[a hh 5 KS S 3 xy(S 4 + a 2 bb z 2 )[-2a 2 bh A KS S + v{-A KS (3 + up)} + S(a 8 bb A KS z 4 
+at h A KS z 2 (2S 4 + 5 KS x 2 z 2 ) + a 4 bb A KS S 4 [S 4 + 25 KS (x - y)(x + y)z 2 \ 
+S 7 y 2 u[- A KS P(1 + C) + M + a 2 bh S 3 [-2A KS S 5 y 2 (l + ()+ S KS (A KS S 5 x 2 
-A KS y 2 z 2 f3v + y 2 z 2 u 2 p)} + A KS S{2a 3 bh 5 KS xyz 2 (S 4 + a 2 hh z 2 ) + [2a 2 hh S 8 + 2a e hh z 4 
+a 4 bh (AS 4 z 2 + 8 KS y 2 z 4 ) - S 8 y 2 (l + ()]S + (S 4 + a 2 bh z 2 ) 2 S 3 })}} , (57) 



m = A 2 KS z^ 2 [ S4 ^ + S ^ AKS ~ + S 4 (at + 5 2 ), C P 2 

+S 3 (A^ 5 /3 - ^)}{,S 3 y(-2a 2 h AKsS + v(-A KS fi + up)) 
+2a bb A KS S[x(S 4 + a 2 bh z 2 ) + a bh ^ 2 ,S]})], (58) 
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533 = is 2 <S> 2 {z\A 2 KS S 6 l4S 2 (l + C) +45/3(1 + () + f3 2 (l + ( - S KS )} 
-2a 6 hh S[2(-l + 5 K s)S 2 (-A KS + S 2 ) + 25 KS A KS z 2 

+ (-1 + 6 KS )S(-A KS + S 2 )f3]p - a 8 hh (-l + 5 KS )S 2 p 2 + 2a 2 hh A KS S 3 {A K s[25 KS z 2 (3 
-2S 2 p(-2 + S KS - 2C) + 4S 3 (1 + C) + SP 2 (1 - S KS + Q] + S\5 KS f3 - p)p} 
+at{A 2 KS {-A5 KS z i + 45 KS Sz 2 p + 4S 3 /3(1 - S KS + 

+S 4 (4 - 4S KS + 4C) + S 2 (85 KS z 2 + p 2 - 5 KS p 2 + P 2 ()) + AA KS S 3 [(-2 + 5 KS )S 2 
Sksz 2 + (-1 + S KS )Sp]p - (-1 + S KS )S 6 p 2 }) - 2[2A 2 KS (S 4 + a 2 hh z 2 ) 2 
+A KS S(a 2 bh + S 2 )z 2 (2{S 4 + a 2 J-(-l + 6 KS )S 2 + 5 KS z 2 }} - (-1 + S KS )S(a 2 hh 
+S 2 )f3)p + a 2 h (-l + 5 KS )(a 2 hh S + S 3 ) 2 z 2 p 2 }S 2 

-(-1 + 5 KS )S 2 (a 2 hh + S 2 ) 2 z 2 p 2 S 4 }, (59) 



and g a/3 = gp a . Also, 



.2 , 52 



^ = < h + 5 2 , (61) 

p = 2S + /?, (62) 

= /? 2 (^ + ag h cos 2 (g)) 
(/? 2 + 4S 2 )A xs ' 

$ = . (64) 

2A KS {S± + a 2 z 2 )vS 
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P"J qMi- ' 1 ! R=0 .01 R t +0 . 99R st f) qu.^0.45; R=0.0lR t tO.93R ; 



Figure 1: NPV maps for a bh = m bh y / 3/4 for the Kerr-Newman metric. See the text for an 
explanation of the colour mapping, (a) At the surface R = 0.99/2+ + 0.01i?s + with q hh = 0.1; 
(b) at the surface R = 0.5R + + 0.5Rs + with q hh = 0.1; (c) at the surface R = 0.01i? + + 0.99i?s + 
with q hh = 0.1; (d) at the surface R = 0.99i? + + 0.01i?s + with q hh = 0.45; (e) at the surface 
R = 0.5R+ + 0.5R S+ with q bh = 0.45; and (f) at the surface R = 0MR+ +0.99R S+ with q bh = 0.45. 
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r] gtl.^ - 1 ! S=0 .01 S t +0 . 99s st £) q bh = 0.'35, S = □ . 01 S,+ . 3 9S Ef . 



Figure 2: NPV maps for a bh = m b h\/3/4 for the Kerr-Sen metric. See the text for an explanation 
of the colour mapping, (a) At the surface 5 = 0.995+ + 0.015s + with q hh = 0.1; (b) at the surface 
S = 0.55+ + 0.5Ss + with q hh = 0.1; (c) at the surface S = 0.015+ + 0.995 5+ with q hh = 0.1; (d) 
at the surface 5 = 0.995+ + 0.015,s + with q bh = 0.45; (e) at the surface 5 = 0.55+ + 0.55,s + with 
q bh = 0.45; and (f) at the surface 5 = 0.015+ + 0.995 5+ with q hh = 0.45. 
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